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TECH NICAL NOTE NO . 1800 

LAMINAE MTYTN O OF A COMERESSIBIE FLUID 
By Dean E. Chapman 


SUMMARY 

A theoretical investigation of the velocity profiles for laminar 
mixing of a compressible— flui<l stream with a region of fluid at rest 
has been made assuming that the Prandtl n umb er is unity. A method 
which involves only quadratures is presented for calculating the 
velocity profile in the mixing layer for an arbitrary value of the 
free-stream Mach number.' 

Detailed velocity profiles have been calculated for free— stream 
Mach numbers of 0, 1, 2, 3* and 5. For each Mach number, velocity 
profiles are presented for both a linear and a 0.76-power variation 
of viscosity with absolute temperature . The calculations for a linear 
variation are much simpler than those for a 0. 76-power variation. 'It 
is shown that by selecting the constant of proport ionality in the 
linear approximation such that it gives the correct value for the 
viscosity in the high— temperature part of the mixing layer, the 
resulting velocity profiles are in excellent agreement with those 
calculated by a 0.76-power variation. 


INTRODUCTION 

The velocity profile for turbulent mixing at constant pressure 
of an incompressible stream with a dead— air region has been calcu- 
lated by several investigators, principally Tollmien (reference l) . 
These calculations agree well with the available experimental data, 
although the conventional assumptions regarding the mixing length of 
a turbulent flow have since been shown by experiments to be incorrect 
(reference 2). The many difficulties encountered in making precise 
turbulent— mixing calculations are, of course, a consequence of the 
extremely complicated mechanism govern ing all turbulent flows. In 
contradistinction to the case of turbulent mixing , the mechanism 
involved in laminar mixing is known, and the relation between 
stresses and velocity gradients for laminar flow can be written 
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explicitly in. a relatively simple form. The velocity profiles for 
laminar mixing, however, apparently have not as yet been calculated 
even for the case 'of incompressible flow. It is the purpose of the 
present paper to calculate the velocity profiles for laminar mixing 
starting with zero boundary— layer thickness of an air stream of 
arbitrary temperature and a dead— air region also of arbitrary temper- 
ature. In cases where a laminar boundary layer of appreciable 
thickness exists at the point where ml xlng begins, the results given 
herein are not directly applicable in the initial part of the mixing 
region. For such cases, it is necessary to make some supplementary 
approximation In order to apply the results. 

Since the practical applications of laminar-mixing phenomenon 
usually Involve the flow of a compressible fluid, the present analy- 
sis includes the effects of compressibility. Examples of typical 
flows wherein laminar mixing occurs can be found in the flow of 
small-scale jet pumps, in the flow behind the intersection of shock 
waves of unequal strength, and in the flow immediately behind the 
base of a body which has a laminar boundary layer. 1 


SIMBOIS AND ROTATIONS 

C constant of proportionality between viscosity and temperature 

Cp specific heat at constant pressure 

k coefficient of. heat conduction 

L characteristic length 

M Mach number 

p static pressure 

Pr Prandtl number 



1 The present analysis was undertaken as part of an investigation of 
this latter problem, and originally appeared as Appendix B of a 
thesis "Base Pressure at Supersonic Velocities," submitted to the 
California Institute of Technology, June 19^8. The results of 
some supplementary computations not given in the thesis have been 
added for sake of completeness in the present report. 
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S Sutherlands constant, 1 approximately 216° F for air 

T absolute temperature 

Ujo free— stream velocity 

u,v velocity components in x,y directions, respectively 
x,y Cartesian coordinates 

8 thickness of mixing layer, taken betveen points where the 
velocity is 0,01 and 0.99 of the free— stream velocity 

co exponent of viscosity variation with temperature 

p mass density 

p coefficient of viscosity 

V kinematic coefficient of viscosity 

ijr stream function 

£ dimensionless independent variable 

Subscripts and Superscripts 

* dimensionless variables as defined in equation (l4) 

« free— stream conditions 

o stagnation conditions of the free stream 

d conditions in the dead-air region 


BASIC EQUATIONS AND ASSUMPTIONS 
Basic Equations 

A schematic illustration of the flow under consideration is 
shown in figure 1. In order to make the laminar-mixing process 
amenable to calculation, the usual assumptions are made that the 
layer affected by viscosity is thin and has zero pressure gradient. 
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Under these conditions the formal procedure for estimating the order 
of magnitude of the various terms in the complete Navier-Stokes 
equations for viscous compressible flow can he carried through in 
precisely the same manner as is done in the classical (Prandtl) 
treatment of laminar hound ary-layer flows. In so doing, the dynamic 
equation for the x direction reduces to the familiar boundary- 
layer momentum equation. 

pu si + pT sr ^ (“5; ; (1) 

while the dynamic equation for the y direction reduces to zero on 
both sides. In passing from the Navier— Stokes equation to (l) it 
is to be noted that the usual boundary— layer assumption b/x < < 1 
is violated in the immediate vicinity of point 0 (fig. l) just as in 
the case of boundary-: layer flow over a plate. 

By employing the same considerations on order of magnitudes as 
were used for the complete Navier-Stokes equations, the complete 
differential equation representing the balance of energy in viscous 
compressible flow reduces to 


d( c p T) d( c p T) 

pu — ^ — + pv — vf— 
ox oy 




( 2 ) 


which is, of course , the usual energy equation for laminar boundary- 
layer flow. In addition to equations (l) and (2), the equation 
expressing conservation of mass is needed: 


d(pu) . d(pv) 


(3) 


For a given gas the variation of p and Cp with temperature is 
known; hence, the foregoing system of three partial differential 
equations is completed by the addition of the equation of state for 
a region of constant pressure 


JL. £a 

T co P 


w 
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Assumptions 

In order to solve the above system of equations, the following 
assumptions are made: 


1 . 


°P ~ 


constant 


2. Pr = cpp/k = 1 



C is a constant depending on T^ and T^. 


The second of these assumptions if often made in calculating boundary- 
layer flows when only the velocity profile is desired and not the 
thermal characteristics. The difference between the boundary— layer 
velocity distributions for Pr = 1 and Pr = 0.73 is small, as is 
clearly shown by the numerical results of Emmons and Brainerd 
(reference 3). Since the mixing— layer and boundary— layer flows 
differ only in the boundary conditions and not in the differential 
equations, the effect of assumption 2 may be expected to be similar 
in the two types of viscous flow. At moderate supersonic Mach 
numbers, the use of Pr = 1 for air does not introduce more than 
1— or 2-percent error in the boundary— layer velocity profile; and 
hence, for all practical purposes, the mixing— layer velocity distri- 
bution calculations for Pr = 1 should be sufficiently accurate for 
air. 


Assumption 3 needs some explanation since the introduction of a 
constant C differing from unity in the approximate relation between 
viscosity and temperature apparently has not been used in previous 
work. Usually C is taken as unity, and in such cases the approxi- 
mation p/poo = (t/Toq) 1 ^ gives reliable results for a fixed cd, 
provided the free— stream temperature is restricted to a certain 
range. By introducing the factor C, the approximating equation can 
be made to give the same value as an exact equation at any desired 
temperature in the mixing layer regardless of T ra or to. Assuming 
that Sutherland 1 s equation. 


_p = / T_W S T w + S 
Poo \ Too / T + S 


(5) 


is exact for the variation of viscosity with temperature, then the 
approximate equation can be made exact at any given temperature T 
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"by means of the relation 


C 



/ t _\ 3/2 To ° + S 

\T J T + S 


( 6 ) 


In particular, if the approximating equation is linear in temperature 
(a) = 1, thereby greatly simplifying the boundary— layer equations) and. 
the viscosity is matched, at the temperature T d then the above 
equation gives 


C = 



Tco + S 

T d + S 


(7) 


for the constant C. 


SOLUTION TO BASIC EQUATIONS 

As was first pointed, out by Prandtl in reference b, and. later 
used, to advantage by Busemann and Crocco (references 5 and 6, 
respectively), the consequence of the assumption Pr = 1 when 
applied to boundary— layer flow is that the temperature becomes a 
function only of the velocity. Hence 

CpT = f(u) (8) 

By substituting this relation Into equation (2) and using equation (l) 
in conjunction with the assumption Pr = 1, it follows that the 
energy equation Is automatically satisfied if the function f(u) 
satisfies the ordinary differential equation 

^|+1 = 0 ( 9 ) 

du 

Integrating this equation, using the boundary conditions, 

T = Tto for u = U^, 

T = T d/ for u = 0 (10) 

gives 

f(u) = c p T = CpT d - —■ + j^c p ( T^ - T d ) + 2sd J ( 11 ) 
as the relationship between velocity and temperature . Since the 
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temperature determines the density, equation (ll) also provides a 
means for calculating the density as a function of the velocity. 

Following the method first given for incompressible flow by 
von Mises (reference 7) and later used for compressible flow by 
von Herman and Tsien (reference 8), a transformation is made to a 
new set of independent variables (x,ijr) , where is the stream 
function. By using ijr as one of the independent variables, the 
continuity equation (3) is identically satisfied, and the velocity 
components are given by 


p 3^ 

u = — — , v 
p 3y 


p 3^ 

00 

p 3x 


( 12 ) 


Since the requirements of conservation of energy are fulfilled by 
equation (11), and conservation of mass by equation (12), the only 
equation now remaining to be satisfied is the momentum equation (l). 
If a transformation were made to a completely new set of independent 
variables ( s , ijr) , the transformation formula would be 


3y 3y cfy 3y 3 s p^ cty 3y 3s 

_3 3i{r _3_ + — — pv 3 | 3 b 3 

3x 3x 3ilr 3x 3s p 3i|r 3x 3s 

- 00 


Setting s = x. 


3s , 3s .. 

— = 0 and v— = 1 
oy ox 


so that the transformation formulae are 2 

3 pu 3 

3y p 3i|f 


/V \ 

V3x/ y 


= -91 


3t 



2 The variables held constant in a differentiation process are 

explicitly indicated in those cases where ambiguity could result 
if the subscript notation were not used. 
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It follows that 


and 


pu 



+ pv 




>(tl -KD, 


'Cty'x. Poo '“ T ' X 
Hence the mome ntum equation (l) in the (z,ijf) system becomes 


8u S f p chiN 

p — = — I un — ) 

80 8x \ p 8i|r/ 

00 


(13) 


This can he put in dimensionless form hy introducing the variables 


„ u 

u * -tr 

00 


T 

- J- 
T 

00 


x* = ^ 


M* = Ji = CT* 40 

Poo 


p* = — 


oo 


t* = 


i r 


V^TJooLC 


(HO 


Except for the parameter C appearing in the definition of ^ ijr* 
and p*, these variables are the same as those uBed by Harman and 
Tsien (reference 8) . Remembering that T 0 *> the free-etream total- 
temperature ratio is given by 


T 

T * = — = 
o t 
00 


1 + 


7-1 


M. 


then the relation (equation (ll)) between temperature and velocity 
can be written as 

T* = Td* - M^u* 2 + (T 0 *-Td*)u* (15) 


The momentum equation (13) becomes, using p*T* = 1, 

du* _ d f 

tx* tty* \ tty* J 


( 16 ) 
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This is the basic equation which most be solved. The boundary 
conditions of the problem are such that no boundary layer exists at 
the point where mixing first begins. Under these conditions the 
velocity profiles will be similar at all points downstream of the 
origin, and hence the velocity u* will be a function only of some 
dimensionless variable £ . This dimensionless variable must involve 
both ty* and x*, and must be zero at the origin of coord inat es 
since the mixing— layer thickness is also zero at the origin. 
Therefore, let 


£ = 7* a x* b 

where a and b are pure numbers which must be determined by the 
condition that both sides of equation (l6) for u* are functions 
only of the single variable £. Setting 

gCOmu*!*^ 1 (17) 

then the right side of equation (l6) can be written as 



from which it is obvious that in order for the right side of 
equation (l6) to be a function only of £, it is necessary that 
a = 1. With # = and a = 1, the partial differential 

equation (l 6) reduces to the ordinary differential equation 


b x* 


-eb— l 


f du* d 

c ir = i r 




t = 




*/ UooVooXC 


( 18 ) 


and the ordinary differential equation for the velocity distribution 
now reduces to 


£ du* d / du*\ 

2 "dT = dTV S dT / 


(19) 


Equation (19) is the same differential equation that was 
obtained for boundary-layer flow -in reference 8. It is a nonlinear 
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differential equation since g depends upon the velocity u*. This 
equation, however, can easily he converted into an integral equation 
which can he solved hy the method of successive approximations. The 
conversion is made hy temporarily assuming that g is a known 
function of £ (instead of u*) and formally applying the standard 
methods for solving first— order li n ear differential equations. The 
result is 


u* = Ci 


■f. i« 


where 


-/? k 


'o 2g 


E = e 

The boundary conditions are 

u* = 1 at £ = 


written as 


ii* = Cxjf | d£ + Uo* 


+ C2 

(20) 

d£ 

(21) 

CD 

(22a) 

— CD 

(22h) 

, = 0, equation (20) 

can he 

+ Uq* 

(23) 


The constant Ci must satisfy two requirements 


l-^o* 



= Ci = 



( 24 ) - 


Equation (23) is an integral equation for u*, since hoth E 
and g are functions of u*. By simply estimating a reasonable 
solution for u* as a function of £, a first approximation iu* to 
the true solution can he calculated from 

1U* = Cl P 2ld£ + oUo* 

^o 

The zero-order approximations qF and og can he calculated 
directly from 0 u* hy using equations (15), (17)> and (21). If this 
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process is repeated until a given approximation is the same as the 
preceding one (to the degree of accuracy desired), and equation (24) 
is simultaneously satisfied, then the solution to the problem is 
obtained. The iteration process turns out to he rapidly convergent, 
requiring two or three iterations to obtain the function u*(£) 
accurate to within 1 percent, and about four or five iterations to 
obtain u*(£) accurate to within a- few tenths of 1 percent. 

In order to change the function u*(£) back to the physical 
coordinates (x,y) a simple quadrature is necessary. By definition 
of the stream function, 

»*(*.),*■ 

or 

gfS fi&L dx + J UooVcXC d£ = ^ dy - HI dx 

00 00 


Hence, with x held constant, integration gives 


from which u/u<» 



(25) 

can be determined. It 


is to be noted that no graphical or numerical differentiations are 
needed at any point in the above iteration process, only quadratures 
are required. 


As is evident from equation (l6), the assumption p* = CT*, 
that is, £0 = 1, makes the momentum equation (l6) independent of 
temperature, and hence density. Consequently, with to = 1, the 
solution to equation (l6) in (x,t|c) coordinates is independent of 
Mach number. For zero Mach number, T* = 1, C = 1, and 



u* 


Using this relation the solution in physical (x,y) coordinates is 
obtained from the solution for zero Mach number by substituting 
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equation (15) into equation ( 25 ) to give 





u*a£+ (T 0 *-Td*H 


( 26 ) 


where the integral in the second term on the right side is carried 
out for u*(£) corresponding to M w = 0. 


RESULTS 

Numerical calculations of the velocity distribution have been 
made for the following cases: 

1. v l * = T*° ‘ 76 j Moo = 0, 1, 2, 3, and 5 

2. p* = CT*; Mao = 0, 1, 2, 3, and 5 

The various solutions for case 2 are obtained directly 
(equation (26)) from the solution for M^, = 0 without carrying out 
the laborious iteration process that is necessary to obtain solutions 
for case 1. All numerical results in physical coordinates 

(y VDW r v®x as independent variable) have been calculated for the 
case Ta = T 0 . If the dead-air temperature is radically different 

from the free— stream stagnation temperature the proper velocity 
distribution curves can be obtained by carrying out the integration 
indicated in equation (25 ) } since the function u*(£) in (x,ijr) 
coordinates is independent of the thermal boundary conditions of the 
problem. 

Curves of u*( £) are shown in figure 2 for various Mach 
numbers. The corresponding curves in the physical plane are shown 
in figure 3 for the case p* = T*°* 76 , and in figure 4 for the case 
p* = CT*. In the latter two figures the familiar Blasius curve for 
the incompressible l ami nar boundary— layer flow is shown for purposes 
of comparison. The constant C that is used in figure 4 is 
determined by matc hing the viscosity coefficient at the temperature 
T^* = T 0 * accord ing to equation (7) • 

The particular curves shown in figure 4 apply for Too — 4oo° B. 
Curves for any other temperature level T,,, differ only in the 
constant factor C. 
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CONCLUDING REMARKS 

A comparison is shown in figure 5 which illustrates the good 
agreement between velocity distributions calculated for the two 
approximations, p* = T*° * 76 , and p* = CT*. For general use the 
linear approximation is recommended since it gives results which 
are practically as accurate as the former, yet does not require a 
laborious iteration solution to be worked out for each Mach number - 


In general, the’ -laminar-mixing layer is several times thicker 
than the laminar boundary layer, as is illustrated in figures 3 and 
4 where, for purposes of comparison, the Blasius profile is alBo 
shown. The rate of growth of mixing— layer thickness with increasing 
Mach number is somewhat larger than the corresponding rate of growth 
for a laminar boundary layer. The curves in figure 3 indicate a 
value of roughly 

®k 


^MfO 


= 1 + 0.11 Meo 2 


for rate of growth of the mixing layer; whereas the corresponding 
value for a laminar boundary layer (see reference 9, for example ) is 
approximately 

8 M p 

K " 1 * °- 09 M " 

This difference is to be expected since a larger percentage of low 
density air exists in a mixing layer than in a boundary layer. 

The foregoing statements, which indicate an increase in mixing- 
layer thickness with increasing Mach number, are based on the assump- 
tion that the Reynolds number (XJ^k/v^) is held constant while the 
Mach number is varied. In most experimental apparatus the Reynolds 
number changes considerably with a variation in Mach number. 
Consequently, depending upon the particular experimental method 
employed, the observed rate of mixing in the x direction may be 
either increased or decreased if the Mach number is increased. 
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Figure 4. -Velocity distribution for linear variation of viscosity 
with temperature. 





